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ABSTRACT 


An  all  digital  simulation  of  the  existing  boiler  and  control  system 
of  CVA-63  was  carried  out.   A  significant  savings  in  computer  time  over 
previous  simulations  was  realized. 

The  parameter-plane  method  was  used  to  search  for  new  operating 
points  for  each  of  the  system  minor  control  loops.   Individual  loop 
transient  response  was  tested  with  a  simulation  program  to  demonstrate 
improved  response  with  new  operating  points. 

The  entire  system  response  was  then  investigated  to  demonstrate  the 
overall  effects.   Recommendations  were  made  for  new  controller  settings 
along  with  component  variations  in  order  to  improve  over-all  system 
response. 
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I.   INTRODUCTION 

Naval  boilers  differ  from  most  other  stationary  type  boilers  in 
that  they  are  called  upon  to  answer  widely  varying  load  demands  in  a 
minimum  amount  of  time.   Feedwater,  fuel  oil  and  combustion  air  must  be 
supplied  in  varying,  but  accurate,  amounts  to  maintain  sufficient  steam 
flow  as  required  by  ship  operation.   There  are  maneuvering  conditions 
which  might  require  boiler  load  to  change  from  ten  or  twenty  percent  to 
over  one  hundred  percent  in  a  matter  of  a  few  seconds. 

Automatic  combustion  control  systems  have  been  developed  by  various 
manufacturers  to  control  the  three  inputs;  oil,  air  and  water.   Most  of 
the  ACC  systems  up  until  the  present  have  been  of  the  pneumatic  variety; 
i.e.,  pneumatic  sensing  and  valve  control.   Individual  components  of 
systems  of  this  type  range  from  simple  bellows  type  sensors  to  the  ex- 
tremely complex  masses  of  springs,  levers,  valves,  gages  and  nozzles 
found  in  pressure  relays  and  transmitters. 

This  thesis  will  investigate  the  automatic  combustion  control  system 
including  a  main  boiler  as  found  on  USS  KITTY  HAWK  (CVA-63).  An  all- 
digital  computer  model  will  be  developed  for  testing  and  analysis  purposes. 
Each  of  the  three  major  control  loops,  fuel  flow,  air  flow  and  feedwater 
flow,  will  be  investigated  independently  and  the  effects  of  any  changes 
or  improvements  will  be  tested  in  the  system  simulation.   No  radical 
changes  to  the  overall  system  will  be  proposed  but  readjustment  of  vari- 
able parameters  within  the  existing  system  will  be  stressed  with  emphasis 
on  improved  overall  system  performance. 

A.   THE  BOILER 

The  boiler   to  be  modeled   is   a  Foster-Wheeler  D-type  as    found   on 
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CVA-63.   The  boiler  consists  of  an  upper  steam  drum  connected  by  tubes  to 
the  lower  water  drum,  screen  wall  and  waterwall  headers.  A  vertical  super- 
heater is  integral  to  the  boiler.   This  unit  is  capable  of  generating 
318,000  lb./hr.  of  steam  at  120  percent  overload  conditions.   Cruising 
conditions  are  considered  to  be  117,500  lb./hr.  of  steam  which  is  44.7 
percent  of  full  boiler  load.   See  Appendix  A  and  Ref.  1  for  more  design 
and  performance  data.   Two  boilers  are  operated  together  and  controlled 
by  one  ACC  unit  due  to  physical  location  on  shipboard.   Boiler  transfer 
functions  were  obtained  from  Ref.  2  and  were  valid  for  cruising  condi- 
tions only.  For  this  reason,  small  variations  in  boiler  load  were  used 
and  the  model  developed  was  of  the  perturbation  variety.  When  data  on 
other  conditions  of  load  become  available,  the  model  may  be  easily  modi- 
fied to  suit  the  design  engineer. 

B.   AUTOMATIC  COMBUSTION  CONTROL 

The  automatic  combustion  control  system  encompasses  two  classifica- 
tions of  pneumatic  control,  combustion  control  and  feedwater  control. 
While  the  feedwater  portion  is  not  essential  for  combustion,  it  is  usually 
included  in  any  discussion  of  ACC  systems  since  it  is  very  necessary  for 
operation  of  the  overall  system. 

The  primary  function  of  the  combustion  control  portion  of  the  system 
is  to  maintain  constant  boiler  superheater  outlet  pressure  under  all  load 
conditions.   It  must  also  insure  that  fuel  and  air  are  mixed  in  the 
proper  proportions  to  maintain  optimum  combustion  efficiency.   Normally, 
a  safety  feature  is  p  ovided  with  the  system  whereby,  in  the  event  of  loss 
of  control  air  pressure,  all  valves  return  to  a  neutral  position  which 
reduces  boiler  firing  to  a  stable  minimum. 
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The  feedwater  control  portion  of  the  system  must  maintain  a  balance 
between  feedwater  flow  to  the  boiler  and  steam  flow  from  the  superheater. 
It  must  also  maintain  the  water  level  in  the  steam  drum  at  a  normal 
steaming  level  under  all  conditions  of  boiler  load.   This  level  control 
is  a  most  important  function  since  high  water  level  can  mean  possible 
liquid-state  carryover  into  rotating  turbines  and  subsequent  destruction 
of  blading  while  low  water  can  expose  the  ends  of  tubes  within  the  boiler 
leading  to  sagging  and  weakening  due  to  extreme  heat  without  adequate 
circulation  cooling.   Level  control  is  complicated  by  the  phenomena  of 
shrink  and  swell.   When  boiler  load  is  increased  rapidly,  large  amounts 
of  steam  are  removed  from  the  boiler  and  the  water  in  the  steam  drum  will 
surge  to  a  high  level.   The  reverse  is  true  when  the  load  is  suddenly 
decreased.   The  thermodynamic  explanation  of  shrink  and  swell  is  to  be 
found  in  Ref.  3  while  more  discussion  is  found  in  this  paper  in  Section 
II.A.6. 

C.    OPERATION  OF  AUTOMATIC  COMBUSTION  CONTROL  SYSTEM 

The  steam  pressure  transmitter  (Fig.  Id)  must  transmit  a  linear 
pneumatic  signal  directly  proportional  to  pressure  variations  in  the  two- 
boiler  main  steam  header.  At  the  normal  operating  pressure,  1200  psig., 
the  output  of  the  steam  pressure  transmitter  is  9  psig.   The  steam  pres- 
sure controller  is  of  the  set-reset  variety,  also  known  as  proportional 
plus  integral  control.  When  the  steam  pressure  signal  from  the  steam 
pressure  transmitter  decreases,  the  output  of  the  steam  pressure  controller 
increases  proportionally  and  conversely.   The  proportional  band  setting 
of  the  controller  determines  magnitude  of  output  variation  for  a  given  in- 
put and  the  reset  needle  valve  of  the  controller  determines  the  rate 
response. 
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The  output  of  the  steam  pressure  controller  goes  both  to  the  air 
flow  controller  and  to  the  low-pressure  signal  selector.   This  output 
becomes  the  demand  signal  for  the  air  flow  control  loop.   The  air  flow 
controller  compares  this  signal  to  the  actual  air  flow  signal  which  is 
proportional  to  the  flow  of  combustion  air  across  the  boiler  air  regis- 
ters.  The  air  flow  controller  is  also  of  the  proportional  plus  integral 
type  and  its  output  signal  is  sent  to  the  airflow  ratio  relay  and  damper 
actuators.  There  are  two  factors  which  affect  the  quantity  of  air  flow, 
damper  position  and  forced  draft  blower  turbine  speed.   The  ratio  relay 
preceeding  each  blower  controller  insures  parallel  operation  of  blowers 
for  there  can  be  one,  two  or  three  blowers  in  operation  for  each  boiler 
at  any  one  time.   It  has  been  found  that  the  use  of  dampers  is  not  only 
unnecessary  but  that  they  have  actually  contributed  to  some  air  flow  loop 
instability  problems.   Dampers  were  completely  eliminated  on  USS  AMERICA 
(CVA-66)  (See  Ref.  4). 

The  air  flow  transmitter  senses  pressure  differential  between  wind- 
box  and  furnace  and  transmits  a  signal  proportional  to  this  difference. 
This  signal,  representing  actual  airflow  across  the  burner  registers,  is 
sent  to  a  ratio  relay  where  the  important  air-fuel  ratio  may  be  adjusted 
for  optimum  combustion  efficiency. 

The  output  of  the  ratio  relay  goes  to  the  air  flow  controller  already 
discussed  and  to  the  low-pressure  signal  selector.   This  signal  selector 
transmits  the  lower  of  its  two  input  signals,  boiler  load  demand  from 
the  steam  pressure  controller  and  air  flow  signal  from  the  fuel-air  ratio 
relay.   The  lower  of  these  two  signals  is  used  as  a  demand  signal  to  the 
oil  flow  portion  of  the  system  so  that,  under  increasing  load,  oil  sup- 
plied to  the  furnace  can  never  exceed  the  required  air  flow  for  proper 
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combustion.   This  prevents  the  formation  of  black  smoke  in  the  furnace 
when  the  load  is  suddenly  increased. 

The  output  of  the  signal  s  lector  goes  to  a  characterizing  relay 
which  converts  it  to  a  corresponding  non-linear  signal  conforming  to  the 
burner  curve  of  return  oil  pressure  vs.  burner  oil  rate.   Fuel  oil  is 
supplied  at  a  constant  supply  pressure  and  the  amount  burned  in  the 
furnace  is  controlled  by  closing  down  on  the  return  valve  in  the  recircu- 
lating system. 

The  feedwater  portion  of  the  control  system  consists  of  a  thermo- 
static drum  level  sensor,  a  combining  relay  and  the  feedwater  regulating 
valve. 

The  level  sensor  consists  of  a  stainless  steel  expansion  tube  in- 
clined across  the  steam  drum  with  its  center  at  normal  steaming  water 
level.   The  upper  tube  end  is  rigidly  fixed  and  the  lower  end  tube  move- 
ment is  amplified  through  a  lever  system  attached  to  a  pilot  valve  stem. 
Level  sensing  motion  comes  from  contraction  due  to  temperature  difference 
between  saturated  steam  temperature  (upper  half)  and  boiler  water  tempera- 
ture (lower  half).   The  pilot  valve  converts  the  mechanical  motion  into  a 
proportional  air  output  signal. 

A  boiler  steam  flow  signal  and  water  flow  signal  are  combined  with 
this  level  signal  in  the  combining  relay  using  a  combination  of  diaphragms 
and  levers.   The  resultant  signal  is  a  total  water  demand  signal  and 
controls  the  feedwater  regulating  valve.   The  feedwater  valve  positioner 
controls  the  flow  of  feedwater  by  converting  the  demand  signal  to  mechani- 
cal valve  motion. 

Now  that  a  general  system  description  has  been  given,  consider  the 
effects  of  changing  the  boiler  load. 
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1.  System  Action  During  Increasing  Boiler  Load 

The  first  noticable  effect  on  increasing  boiler  load  is  a  drop 
in  superheater  outlet  pressure.   This  is  reflected  by  a  decrease  in  the 
output  of  the  steam  pressure  transmitter  which  causes  a  corresponding  in- 
crease in  the  output  of  the  steam  pressure  controller. 

This  represents  an  increasing  air  demand  signal  which  is  greater 
than  the  actual  airflow  signal.   This  causes  an  increase  in  blower  speed 
and  opens  dampers  further  until  flow  matches  demand. 

When  load  increases,  the  signal  transmitted  by  the  low-pressure 
signal  selector  is  the  srriller  air  flow  signal.  As  the  air  flow  return 
system  receives  an  increasing  signal,  the  fuel  oil  return  valve  is  closed 
proportionately  causing  more  oil  to  be  burned  in  the  furnace.   This,  in 
turn,  causes  an  increase  in  steam  pressure  until  equilibrium  is  once  again 
reached. 

The  effects  of  boiler  water  swell  are  reduced  by  the  action  of 
the  combining  relay.   As  load  increases,  the  steam  flow  signal  increases 
to  a  new,  higher  fixed  position.  The  swelling  water  level  produces  a 
false  signal  (actually,  more  water  is  needed,  not  less  as  indicated  at 
first)  which  is  balanced  by  the  steam  flow  signal.   Relay  settings  deter- 
mine sensitivity.  As  the  level  swell  boils  away,  the  level  signal  will 
decrease  to  its  true  value  and  a  demand  will  be  made  for  more  water.  As 
water  flow  increases  to  a  higher  value,  its  signal  will  offset  the  steam 
flow  signal  and  equilibrium  will  be  reached  once  again. 

2.  System  Action  During  Decreasing  Boiler  Load 

When  boiler  load  is  decreased,  the  superheater  outlet  pressure 
increases  causing  an  increase  in  steam  pressure  signal  from  the  steam 
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pressure  transmitter.   This  increasing  signal  causes  a  proportional  de- 
crease in  steam  pressure  controller  output.   Since  air  flow  is  now  greater 
than  necessary,  the  output  of  the  minimum  signal  selector  becomes  the  lower 
steam  pressure  signal,  decreasing  fuel  oil  demand.   The  air  flow  response 
lags  behind  at  this  point,  providing  an  excess  of  combustion  air  until 
blowers  slow  and  dampers  close. 

The  decrease  in  fuel  oil  demand  opens  the  return  valve  and  the 
pressure  returns  to  equilibrium  due  to  the  smaller  quantity  of  fuel  oil 
being  burned  in  the  furnace. 

The  water  flow  portion  will  act  oppositely  to  the  manner  in 
which  it  acted  for  a  load  increase.   Decreased  steam  flow  demand  will  be 
partially  (or  fully,  depending  on  adjustment  of  combining  relay)  offset 
by  shrink  in  boiler  water  level  signal.  As  this  shrink  effect  dies  out 
and  the  drum  level  starts  to  increase,  the  output  of  the  combining  relay 
will  increase  and  the  feedwater  regulating  valve  will  close.  As  water 
flow  decreases,  the  water  flow  signal  offsets  the  steam  flow  signal  in 
the  combining  relay  and  equilibrium  is  reached. 

D.    IMPROVEMENTS  NEEDED 

Control  of  a  boiler  and  its  associated  auxiliary  equipment  is  compli- 
cated considerably  by  the  interaction  present  among  its  control  variables. 
The  representation  of  the  boiler  and  control  system  in  Figs.  1(a)  through 
1(d)  and  the  interconnection  diagram  of  Fig.  2  demonstrate  the  coupling 
which  exists.   All  three  major  system  variables;  steam  flow,  oil  flow, 
and  water  flow  affect  both  boiler  steam  pressure  and  boiler  water  level 
in  varying  ways.   The  transfer  functions  seem  complicated  enough  at  first 
inspection  but  when  one  stops  to  consider  that  these  are  valid  for  cruising 
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conditions  only  and  that  at  varying  boiler  loads  they  become  dynamic 
functions  of  the  load  and  rate  of  load  change  the  control  problem  becomes 
overwhelming. 

Ideally,  if  the  interaction  between  inputs  and  outputs  could  be  re- 
moved or  balanced  in  some  way,  the  control  problem  would  be  greatly 
reduced  in  scope.   In  the  case  of  the  boiler,  physically  forcing  non- 
interaction is  impossible.   The  alternative  then,  is  to  design  a  non- 
interacting  controller.   This  has,  in  fact,  been  done  on  paper  (See  Ref. 
5),  but  the  controller  transfer  functions  which  result  for  this  system 
are  not  physically  realizable  and  could  not  be  simulated  using  the  digital 
simulation  language  to  be  discussed  in  the  next  section.   See  Appendix  B 
for  design  information  on  the  non- interacting  controller. 

Any  improvement  then,  in  control  of  this  boiler,  is  limited  to 
complete  redesign  of  the  entire  system,  or  optimization  of  the  existing 
system. 

Alternative  control  systems  have  been  proposed.   Digital  control 
(Refs.  6  and  7)  and  analog  fluidic  control  (Ref.  8)  are  two  of  the  better 
known  alternatives.   The  first  was  investigated  by  Naval  Electronics  Labor- 
atory Center  but  this  investigation  was  terminated  in  June  1968  due  to 
unfavorable  reaction  by  Naval  Ship  Engineering  Center.  Analog  fluidic 
control  is  under  test  and  development  presently  at  NAVSEC,  Philadelphia. 
It  is  safe  to  assume  that  this  testing  and  design  will  come  to  fruition 
in  a  new,  fast,  state-of-the-art  control  system  for  ships  of  the  future. 

What  are  we  to  do  with  the  ships  of  today?  Economics  forbids  the 
complete  replacement  of  several  hundreds  of  ACC  systems  on  ships  which 
have  a  life  expectancy  of  between  twenty  and  thirty  years.   Hence,  let  us 
look  into  the  possibility  of  improving  present  systems  at  a  minimum  of 
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expense.   At  the  very  least,  we  can  hope  to  explore  methods  which  will 
enable  the  future  ship  control  designer  to  design  better  systems. 

Exactly  what  is  needed  in  the  way  of  improvements?  Basically,  it 
is  a  faster  responding  system  with  less  overshoot  of  variable  response. 
For  example,  Ref.  9  states  that  specifications  on  several  ships  required 
boiler  water  level  to  remain  within  +  4  inches  of  set  point  during  linear 
changes  in  boiler  load  from  ten  percent  to  eighty  percent  in  45  seconds. 
Subsequent  shipboard  operation  produced  high  water  conditions  out  of 
sight  in  the  boiler  gage  glass.   This  was  due  to  maneuvering  the  ship  with- 
in the  load  specifications  but  in  only  10  seconds  instead  of  the  specified 
45  seconds.  Although  the  system  would  probably  meet  the  original  specifi- 
cations, response  could  probably  be  improved  somewhat  through  a  judicious 
"tuning"  of  the  control  system. 

This  thesis  investigates  the  CVA-63  system  and  makes  recommendations 
for  readjustment  of  controller  settings  or  replacement  of  component  parts 
as  necessary  for  better  (faster)  overall  system  response.   The  parameter 
plane  method  of  optimization  is  used  and  it  is  hoped  that  even  if  system 
response  is  not  dramatically  improved,  the  methods  used  will  prove  en- 
lightening for  others  who  may  follow. 
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II.   THE  SIMULATION 

The  boiler  and  automatic  combustion  control  system  were  simulated 
utilizing  the  IBM  developed  Digital  Simulation  Language  (DSL).   DSL  is 
a  non-procedural,  problem-oriented  language  which  will  accept  problems 
in  ordinary  differential  equation  form  or  as  analog  block  diagrams.   Since 
an  excellent  block  diagram  of  the  CVA-63  system  was  available  (Ref.  2), 
this  program  provided  a  relatively  simple,  efficient  means  of  simulation. 
References  10  and  11  describe  the  DSL  program  in  more  detail.   The  vari- 
ables indicated  in  parenthesis  refer  to  the  system  diagrams,  Figs.  1(a)  - 
1(d)  and  the  interconnection  diagram,  Fig.  2. 

A.   BOILER 

The  seven  transfer  functions  appearing  in  Fig.  1(a)  are  the  result 
of  tests  conducted  at  NBTL  (Ref.  2).   These  functions  are  the  result  of 
a  combination  of  both  pulse  testing  (Ref.  12)  and  of  periodic  input  fre- 
quency response  analysis. 

1.  Steam  Pressure  /  Steam  Flow  (PDS/GS) 

This  response  essentially  represents  a  mass  balance  situation 
and  is  accurately  represented  by  a  straight  integration  at  cruising  condi- 
tions. 

2.  Steam  Pressure  /  Fuel  Oil  Flow  (PDF/GF) 

This  can  be  represented  by  a  first-order  lag,  first-order  lead, 
and  a  straight  integration,  all  series  connected. 

3.  Steam  Pressure  /  Water  Flow  (POW/GW) 

A  second-order  lag  (complex  poles)  and  series  connected  integra- 
tion represent  the  effects  of  feedwater  flow  on  steam  pressure  at  cruising 
conditions. 
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4.  Water  Level  /  Water  Flow  (LW/GW) 

As  in  1.  on  the  preceeding  page,  a  straight  integration  ac- 
curately represents  this  simple  mass  balance  situation. 

5.  Water  Level  /  Oil  Flow  (LF/GF) 

Pulse  tests  by  NBTL  yielded  a  transfer  function  consisting  of 
integration  and  series-connected  lag-lead. 

6.  Water  Level  /  Steam  Flow  (LST/GS) 

The  response  of  water  level  to  steam  flow  rate  is  shown  by  the 
parallel  combination  with  outputs  LSA  and  LSB  in  Fig.  1(a).   Water  level 
in  the  steam  drum  should  integrate  proportionally  to  the  mass  balance  of 
steam  flow-water  flow  difference.   This  integration,  however,  is  of  low 
natural  frequency  and  is  initially  overshadowed  by  the  shrink/swell 
phenomenon  discussed  previously.  Fig.  3  gives  the  response,  LST,  to  a 
step  input.   Fig.  4  give"  the  same  response  to  the  triangular  pulse  input 
used  for  response  tests  in  the  system  (Fig.  5).   It  can  be  seen  that  the 
effects  of  the  integration  are  not  present  for  the  first  several  seconds. 

B.    FUEL  OIL  FLOW  CONTROL 

Also  called  steam  pressure  control  loop,  the  fuel  oil  flow  control 
loop  is  a  closed  control  loop  consisting  of  the  following  elements; 

1.  Steam  Pressure  Controller  (PZP/PM) 

This  is  a  typical  proportional  plus  integral  controller.  Values 
for  proportional  and  integral  gain  were  measured  by  NBTL  using  sinusoidal 
testing  techniques  and  are  shown  in  Fig.  1(b).   These  gains  are  adjust- 
able and  will  be  utilized  to  set  the  loop  to  secure  "optimum"  performance. 

2.  Low  Pressure  Signal  Selector  (PM,  PQR/PSF) 

This  acts  as  a  switch  only  and,  as  such,  has  no  measurable  trans- 
fer function.   It  is  conveniently  simulated  in  DSL  by  using  a  Function 

Switch  (see  computer  output). 
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3.  Characterizing  Relay  (PSF/PRF) 

This  has  a  two-pole  response  as  determined  by  NBTL  using  the 
periodic  input  method. 

4.  Fuel  Oil  Control  Valve  (PRF/PF) 

A  quadratic  lag  function  was  obtained  for  this  valve.   A  small 
time  delay  (0.023  seconds)  is  also  present  due  to  valve-stem  stiction. 
This  delay  appears  in  all  of  the  pneumatic/mechanical  valves  and  is  present 
due  to  the  physical  limitations  of  the  valves  themselves,  and  their  in- 
ability to  respond  instantaneously  to  a  demand. 

5.  Steam  Pressure  Transmitter  (POUT/PP) 

Analysis  of  transient  response  tests  by  NBTL  yielded  a  quadratic 
lag  transfer  function  as  shown  in  Fig.  1(b). 

C.   WATER  LEVEL/FLOW  CONTROL 

The  water  control  loop  is  a  "three-element"  system,  responding  to 
changes  in  water  level,  water  flow,  and  steam  flow.   Individual  elements 
(see  Fig.  1(c))  are  as  follows: 

1.  Feedwater  Control  Valve  (PDW/GW) 

The  transfer  function  for  this  valve  was  found  by  NBTL  to  be 
third  order,  consisting  of  a  first-order  lag,  attributed  to  the  capaci- 
tance of  the  pneumatic  diaphragm  chamber,  and  a  second-order  lag  due  to 
the  valve  positioner.  A  0.45-second  pure  time  delay  was  thought  to  be 
due  to  valve-stem  friction  in  the  packing  gland. 

2.  Combining  Relay  (PSA,PWA/PSW) 

Dynamic  measurement  could  not  be  made  by  NBTL  and  this  com- 
ponent is  approximated  by  pure  gains. 

3.  Thermostatic  Drum  Level  Indicator  (L/PL) 

At  cruising  condition  loads,  the  thermostat  may  be  represented 
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by  a  pure  gain.   NBTL  was  unable  to  make  dynamic  measurements  due  to 
limitations  in  the  generation  of  a  suitable  water  level  input  signal. 

D.  AIR  FLOW  CONTROL 

The  air  flow  loop  is  a  typical  feedback  control  system  and  the 
individual  component  transfer  functions  were  derived  by  NBTL  and  are  as 
follows: 

1.  Controller  (POA/PEQ) 

This  controller  is  a  typical  set-reset  controller  with  adjust- 
able integral  and  proportional  gains. 

2.  Ratio  Relay  (PG/POA  and  PQR/PQA) 

The  transfer  function  obtained  by  NBTL  is  as  shown  in  Fig.  1(d). 

3 .  Forced  Draft;  Blower  Actuator  and  Valve  (GB/PQ) 

The  actuator  and  valve  exhibit  characteristic  underdamped 
response  and  time  delay  which  is  normal  with  components  of  this  type. 
Valve-stem  friction  is  not  as  pronounced  as  in  the  feedwater  control  valve 
due  to  decreased  valve-stem  travel  and  the  fact  that  the  medium  being 
controlled  is  steam  instead  of  water. 

4.  Main  Forced  Draft  Blowers  (QA/GB) 

Direct  frequency  response  testing  was  used  by  NBTL  to  derive 
the  transfer  function  shown  in  Fig.  1(d). 

5.  Air  Flow  Transmitter  (PQA/QA) 

This  component  is  represented  by  a  quadratic  transfer  function 
and  appears  in  the  feedback  system  of  the  air  flow  control  loop. 

E.  INITIAL  SYSTEM  TESTS 

Individual  components  were  simulated  and  interconnected  in  DSL  (see 
program  listing)  to  produce  a  system  simulation  for  transient  response 
testing.   Since  transient  response  data  for  this  boiler  and  combustion 
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control  system  are  almost  non-existent,  the  author  had  to  rely  on  personal 
experience  as  Boiler  Officer  on  USS  INDEPENDENCE  (CVA-62)  from  1964  to 
1966.   General  shape  and  accuracy  of  system  response  were  also  verified 
by  J.  W.  Banham,  Jr. ,  head  of  the  Machinery  Automation  Systems  Depart- 
ment of  the  Naval  Ship  Engineering  Center,  Philadelphia  Division. 

Inspection  of  the  completed  simulation  program  on  page  96  will  re- 
veal the  manner  in  which  individual  components  and  loops  are  interconnected, 

It  was  found  necessary  to  put  limitations  on  certain  pneumatic  signals 
to  simulate  the  loading  effects  of  a  controller  or  relay  being  overdriven. 
This  portion  of  the  simulation  appears  at  the  end. 

A  triangular  shaped  input  was  selected  for  transient  response  testing. 
Since  the  model  is  of  the  perturbation  variety,  the  use  of  a  step  input 
was  deemed  impractical  and  unrealistic.   The  triangular  shaped  input  seemed 
to  be  the  closest  meaningful  input  signal  which  most  nearly  resembled  a 
true  impulse.   The  input  is  shown  in  Fig.  5. 

System  responses  to  this  input  are  shown  in  Fig.  6  through  Fig.  10. 
It  must  be  remembered  that  variations  are  measured  from  a  steady  cruising 
condition  value;  i.e.,  a  steam  pressure  variation  of  -8  psig  is  actually 
1192  psig  drum  pressure,  8  psig  below  the  nominal  1200  psig. 

The  value  of  a  component  by  component  simulation  is  readily  apparent 
at  this  stage.  Any  intermediate  signal  may  be  designated  as  output  for 
analysis  purposes.   Components  which  are  overloaded  or  driven  into  satura- 
tion are  easily  discovered  and  individual  components  may  be  replaced  with 
new  or  improved  components  to  test  the  effects  on  the  overall  system 
response. 

The  effect  of  the  underdamped  feedwater  control  valve  on  water  flow 
rate  is  readily  apparent  in  Fig.  6.   Later  transient  response  tests  show 
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the  effects  of  varying  time  delay  in  the  transfer  function  and  of 
changing  one  of  the  gains  in  the  combining  relay. 

The  complete  system  simulation  shows  a  great  time  saving  over  pre- 
vious all-digital  simulations  (Refs.  7  and  13).   Run  time  for  the  DSL 
program  has  a  2:1  correspondence  with  real  or  problem  time.   For  example, 
two  minutes  of  simulated  boiler  operation  takes  one  minute  of  execution 
time  on  the  IBM  360/65.  While  this  depended  somewhat  on  the  size  of  the 
integration  step  and  on  the  number  of  variables  output,  previous  simula- 
tions were  approximately  ten  times  slower  by  comparison. 

Simulation  of  individual  loops  was  used  for  transient  response  testing 
but  the  simple  programs  used  are  subsets  of  the  overall  system  program 
and  do  not  merit  further  discussion. 
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III.   INDIVIDUAL  LOOP  ANALYSIS 

Since  system  complexity  precluded  mass  analysis  of  the  entire  ACC 
system  at  one  time,  it  was  decided  to  separate  the  system  into  three 
minor  loops.   Obvious  choices  for  analysis  were  the  three  major  system 
variables,  combustion  air,  fuel  oil  and  feedwater.  Although  interaction 
is  present  throughout  this  system,  it  was  ignored  for  individual  loop 
analysis  in  the  hope  that  improvements  made  in  each  loop  would  be  signifi- 
cant in  the  whole  system.   This  obviated  the  massive  problem  of  trying  to 
counter  the  effects  of  interaction  to  yield  accurate  loop  responses. 

Several  analysis  methods  exist  for  determination  of  controller  para- 
meters to  permit  "optimum"  performance.   Reference  9  gives  an  excellent 
overview  of  individual  component  relationships  to  system  performance. 
Essentially  however,  due  to  the  demands  of  shipboard  maneuvering,  the 
performance  requirement  can  be  easily  stated: 

The  ACC  system  must  be  capable  of  maintaining  stable  operating 
conditions  under  widely  varying  load  conditions  with  a  given 
set  of  machinery  (boiler,  blowers,  valves,  etc.)  and  it  must 
make  transitions  between  operating  points  in  a  minimum  amount 
of  time. 

This  requirement  is  more  easily  stated  than  met.  The  designer  is  limited 
to  passive  components  and  more  precisely,  to  the  relative  slow  pneumatic- 
mechanical  components  which  are  presently  available. 

Root-Locus,  Bode,  and  Black-Nichols  are  among  analysis  methods 
presently  used  in  the  design  and  analysis  of  systems  of  this  type. 
Reference  14  is  the  first  known  to  apply  the  parameter  plane  technique 
to  the  minor  control  loop  of  a  boiler  ACC  system. 
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The  parameter  plane  is  a  two-parameter  method  which  determines  the 
values  of  two  parameters  (0(andj3)  that  are  required  to  insure  that  the 
characteristic  polynomial  has  roots  at  a  prespecified  location  in  the  s- 
plane.  By  choosing  a  suitable  line  on  the  s-plane  and  mapping  it  through 
certain  algebraic  transformations  (see  Ref.  15)  one  may  obtain  a  curve  on 
the  CX  -p  plane  which  may  be  used  to  select  "optimum"  values  of  0(and  [3 
to  give  a  desired  response.   In  a  dynamic  system  of  this  type  where  speed 
of  response  was  important,  three  mapping  transformations  were  deemed  neces- 
sary.  Constant- CO   lines  ((jt)   is  the  natural  frequency  of  an  approximate 
second-order  system)  map  from  circles  of  constant  radius  on  the  s-plane. 
Constant-  E,   lines  (  ^  is  the  damping  ratio  of  the  approximate  second-order 
system)  map  from  radial  lines  on  the  s-plane.  Constant-  0~  lines  map  from 
real  root  locations  on  the  real  axis  of  the  s-plane. 

One  difficulty  encountered  was  the  mathematical  representation  of  pure 
time  delay  as  present  in  all  three  minor  loops.   It  was  necessary  to  have 
this  delay  in  a  mathematical  form  to  determine  the  characteristic  poly- 
nomial for  each  of  the  loops.   The  Pade  Approximation  of  transport  lag  is 
among  the  most  widely  known  and  was  selected  (Ref.  16).  Varying  degrees 
of  approximations  were  available,  ranging  from  first  to  twelfth  order  with 
accompaning  increases  in  accuracy  with  higher  order.   Since  the  use  of  the 
parameter  plane  was  dependent  on  the  determination  of  the  minor-loop 
characteristic  polynomial,  a  fourth-order  approximation  was  chosen  for  all 
transport  lag  representations.   This  approximation  gave  a  satisfactory 
mathematical  result  and  increased  the  highest  power  of  the  characteristic 
polynomial  by  four.   The  general  form  of  the  approximation  appears  in  Ap- 
pendix C  along  with  the  actual  equation  used  in  each  of  the  loops. 
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Since  hand  calculation  of  mathematical  mapping  transformations 
would  be  prohibitive,  a  computer  program,  developed  at  the  Naval  Post- 
graduate School  under  the  direction  of  Dr.  George  J.  Thaler,  was  used 
for  parameter  plane  graphical  results. 

Algebraic  manipulation  of  transfer  functions  was  also  performed  on 
the  computer.   FORMAC  (FOrmula  MAnipulation  Compiler),  an  algebraic  manipu- 
lation language  (see  Refs.  17  and  18)  was  extremely  useful  for  develop- 
ment of  loop  characteristic  polynomials  and  should  prove  useful  to  the 
control  engineer  faced  with  tedious  algebraic  manipulation. 

All  graphs  presented  in  this  thesis  are  the  result  of  computer  runs 
and  were  drawn  directly  by  a  CalComp  plotter  at  the  computer  facility  at 
the  Naval  Postgraduate  School. 

A.   AIR  FLOW  LOOP  ANALYSIS 

The  air  flow  control  loop  is  a  fairly  typical  feedback  system  and  is 
shown  in  Fig.  11  which  is  a  simplified  version  of  Fig.  1(d).   The  re- 
sulting loop  characteristic  polynomial  for  the  closed  loop  system  is  of 
fifteenth  order  and  is  shown  in  Appendix  C. 

The  parameter  plane  approach  was  used  and  the  selected  variable 
parameters  were  the  proportional  and  integral  gains  of  the  air  flow 
controller.   These  variables  were  thought  to  be  the  most  easily  acces- 
sible and  most  easily  adjusted.   The  proportional  gain  in  the  original 
system  had  been  set  at  1.67  and  the  integral  gain  at  0.024.   For  para- 
meter plane  purposes,  integral  gain  was  designated  beta  ((3  )  and  propor- 
tional gain  alpha  (CX).   The  original  system  operating  point  is  marked 
with  an  X  on  Fig.  12,  the  air  flow  loop  parameter  plane. 
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In  order  to  demonstrate  the  accuracy  and  validity  of  the  parameter 
plane,  several  runs  of  transient  response  were  made  using  a  simple  DSL 
simulation  of  the  system  in  Fig.  11. 

It  was  decided  to  hold  one  parameter  constant  while  varying  the 
other  in  order  to  dramatize  the  effect  of  moving  in  a  straight  line  on 
the  parameter  plane.   In  Fig.  13,  alpha  was  held  constant  at  1.67  while 
beta  was  allowed  to  take  on  values  of  0.01,  0.024,  0.1,  0.2.   Note  that 
the  second  curve;   0(=  1.67,  P =  0.024  is  the  transient  response  of  the 
original  system.   Input  for  all  transient  response  tests  in  this  section 
was  a  unit  step.   Since  motion  along  a  line  of  constant  alpha  constituted 
motion  along  a  line  of  nearly  constant  zeta  (damping),  no  increase  in 
response  time  was  anticipated.   Figure  13  bears  out  this  expectation. 
Motion  along  a  line  of  constant  beta  however,  was  nearly  perpendicular 
to  the  lines  of  constant  zeta,  and  rapid  changes  in  transient  response 
were  expected.   Figure  14  depicts  such  transient  response;  holding  beta 
constant  at  0.024,  alpha  was  allowed  to  take  on  values  of  1.0,  1.67,  2.0, 
2.5,  3.0,  and  3.5.  Again,  for  comparison,  the  01=  1.67,  £  =  0.024 
curve  is  the  response  of  the  original  system.  At  C(=  1.0,  which  is  nearly 
on  the  £>  =  1.0  curve,  the  system  exhibited  the  characteristic  response  of 
a  critically  damped  second-order  system  as  predicted  on  the  parameter 
plane. 

While  no  one  point  could  be  selected  as  the  "optimum1'  operating  point, 
several  qualitative  statements  could  be  made.   It  is  obvious  from  Fig.  13 
that  beta,  at  least  within  the  range  of  investigation,  had  little  effect 
on  the  response  time  of  the  system.   On  the  parameter  plane,  if  alpha  is 
allowed  to  remain  constant,  increasing  beta  will  move  the  operating  point 
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along  a  nearly  constant- c^  path  until  beta  exceeds  about  0.38  at  which 
time  ^  begins  to  decrease.   The  effect  of  variation  of  alpha  was  more 
readily  apparent.   Rise  times  were  considerably  improved  with  alpha  values 
of  2.0  and  2.5  while  overshoot  was  not  considered  to  be  excessive  (Fig. 
14).   On  the  basis  of  this  qualitative  analysis,  the  "improved"  operating 
point  of  0(  =  2.5,   P  =  0.024  was  selected.   This  point  is  marked  with 
a  0  on  the  parameter  plane,  Fig.  12.   Section  IV  demonstrates  the  response 
of  the  entire  system  with  this  new  operating  point. 

B.    OIL  FLOW  LOOP  ANALYSIS 

Since  the  position  of  the  low-pressure  signal  selector  determines 
the  configuration  of  the  oil  flow  control  loop,  a  choice  had  to  be  made. 
During  an  increasing  load,  the  air  flow  system  lags  behind  and  it  is  the 
output  of  this  loop  (PQR  in  Fig.  2)  which  becomes  demand  for  oil  flow. 
Hence,  the  air  flow  loop  is  an  integral  part  of  the  oil  flow  loop  when 
load  is  increasing.   Increasing  load  conditions,  before  the  air  flow  loop 
responds,  the  demand  signal  (PM)  from  the  steam  pressure  controller  is 
lower  and  therefore  becomes  the  demand  signal. 

It  was  decided  to  try  to  improve  the  loop  in  the  simpler  of  the  two 
configurations  since  combination  would  have  led  to  a  characteristic  poly- 
nomial of  the  twenty-eighth  order.   It  was  hoped  that  improvements  made 
to  the  air  flow  loop  would  help  the  response  when  this  loop  was  a  part 
of  the  oil  loop  during  a  load  increase. 

The  simplified  oil  loop  diagram  appears  in  Fig.  15.   The  resulting 
characteristic  polynomial  of  thirteenth  order  was  determined  using  the 
FORMAC  program  and  is  found  in  Appendix  C. 
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As  was  done  with  the  air  flow  loop,  the  parameter  plane  was  plotted. 
This  time,  the  variables  of  the  steam  pressure  controller  were  used  for 
parameters.  Alpha  was  chosen  to  be  the  proportional  gain  and  beta  the 
integral  gain.   The  resulting  graph  is  shown  in  Fig.  16. 

Several  points  on  the  parameter  plane  were  chosen  for  DSL  simulation 
testing.   With  alpha  held  constant  at  6.67,  beta  was  allowed  to  take  on 
the  values  0.001,  0.014,  0.1,  and  1.0.   These  curves  are  shown  in  Fig.  17, 
Note  that   0(=  6.67,  P=»  0.014  is  the  operating  point  of  the  original 
system  as  marked  on  the  parameter  plane  with  an  X.   Three  curves  are  very 
close  together,  the  result  of  moving  vertically  along  a  nearly  constant- 
line  as  before.   The  fourth  curve,  p  =  1.0,  represents  a  point  off  the 
parameter  plane  but  close  inspection  will  reveal  that  the  value  of 
has  been  drastically  reduced  for  this  operating  point.   The  transient 
response  bears  out  this  hypothesis,  demonstrated  by  the  large  overshoot 
of  a  characteristically  underdamped  system.  Again,  with  beta  held  con- 
stant at  0.014,  alpha  was  allowed  to  vary  from  4.0  to  10.0.   These  re- 
sponse curves  are  shown  in  Fig.  18.   The  critical  damping  ratio, ^  =  1.0, 
was  not  approached  with  the  oil  loop  but  extrapolation  of  the  curves  of 
Fig.  16  to  an  alpha  of  2.5  will  lead  to  the  expected  response.   The 
original  system  response  is  shown  for  comparison  as  before  (  #  =  6.67, 
(3  =  0.014). 

Again,  it  was  difficult  to  make  a  quantitative  analysis  to  determine 
an  "optimum"  operating  point.   It  must  be  remembered  that  the  results 
were  valid  only  for  an  increasing  load,  the  loop  being  considerably 
different  during  a  decreasing  load.   A  new  operating  point  was  selected 
on  a  qualitative  basis  at  Ws  8.0,  P  =  0.014  and  this  point  is  marked 
with  a  0  on  the  parameter  plane. 
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C.    WATER  FLOW  LOOP  ANALYSIS 

The  water  flow  control  loop  is  shown  in  simplified  form  in  Fig.  19. 
In  order  to  maintain  true  closed-loop  action  for  the  resulting  analysis, 
the  loop  had  to  be  reduced  as  shown  from  the  loop  in  Fig.  1(c). 

The  integration  in  the  feedback  circuit  made  the  system  character- 
istically type  one.   Hence,  response  to  a  step  input  contained  no  steady- 
state  error,  and  for  unity  input,  unity  output  was  obtained.   As  was  found 
in  the  original  system  response  (see  Fig.  6)  the  response  of  the  water 
loop  left  something  to  be  desired  as  far  as  oscillatory  action  was  con- 
cerned. 

The  loop  characteristic  polynomial  was  determined  using  the  FORMAC 
program  and  is  given  in  Appendix  C. 

Choices  for  parameter  plane  variables  were  somewhat  limited  in  this 
case.   Since  there  was  not  a  proportional  plus  integral  controller  present 
as  in  the  oil  and  air  loops,  other  variables  had  to  be  selected.  A 
simple  inspection  of  the  loop  of  Fig.  19  revealed  two  likely  parameters. 
The  gain  of  the  thermostat  (2.37)  was  selected  as  alpha  and  the  gain  of 
one  of  the  combining  relay  factors  (0.014)  was  selected  as  beta. 

Effects  of  the  i ether  large  (0.45  sec)  time  delay  were  questioned 
in  Ref.  2,  so  several  transient  response  runs  were  made  of  the  entire 
system  and  results  are  shown  and  discussed  in  the  next  section,  IV. D. 

The  parameter  plane  plot  is  shown  in  Fig.  20.   Note  that  the  operating 

point  of  the  original  system,  Of-  -   2.37 ,  p  ■  0.147,  which  is  marked  with 

an  X,  lies  outside  of  the  ^=1.0  curve.   This  indicated  that  the  dominant 

roots  were  not  complex,  but  real.  A  constant-sigma,  or  real  root,  line, 

0"  *  0.52,  passes  through  the  operating  point.   This  indicated  that  one 
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of  the  dominant  roots  lies  at  <T   =  -0.52  on  the  root  locus  plot.   In  order 
to  verify  this  result,  a  computerized  root-locus  plotting  program,  also 
developed  at  the  Naval  Postgraduate  School,  was  used. 

Since  the  root  locus  is  a  one-parameter  method,  it  was  decided  to 
fix  alpha  (thermostat  gain)  at  2.37  and  more  vertically  on  the  parameter 
plane  along  a  straight  line  and  vary  beta  from  0.1  to  1.0.   Since  P=  1.0 
appears  to  be  above  the  C,  =  0.0  curve  on  the  parameter  plane,  it  was  ex- 
pected that  the  roots  would  move  into  the  right  half  of  the  s-plane  on  the 
root  locus,  indicating  instability  of  the  loop.   It  was  further  expected 
that  when  the  constant-alpha  line  crossed  the  S>  =  1.0  line  on  the  para- 
meter plane,  that  the  real  roots  on  the  root  locus  would  come  together 
on  the  real  axis  and  reach  their  breakaway  point.   The  resulting  root- 
locus  plot  is  shown  in  Fig.  21  and  it  was  found  to  fully  verify  the  re- 
sults expected  from  parameter  plane  predictions.  Additional  roots  which 
do  not  appear  on  this  plot  due  to  scale  limitations  were  at  -12.873+j3.8547 
and  -9.3502+J11.811.   System  zeros  were  at  0.0,  12.872+J3.8544  and  9.3502+ 
j 11 .811 .   As  expected,  the  real  roots  which  appeared  at  -1.2176  and  -0.039322 
when  beta  had  the  value  0.1,  moved  together  as  beta  increased.   The  left 
hand  root  was  at  -0.52  at  p  =  0.147  as  predicted  by  the  parameter  plane 
and  both  real  roots  came  together  at  P ■  0.192  (^  =  1.0).   The  roots 
then  became  complex  and  moved  toward  the  right-half  plane,  crossing  over 
at  P  =  0.333  as  predicted  by  the  parameter  plane. 

It  was  expected  that  closed-loop  response  could  be  improved  by  in- 
creasing beta,  as  predicted  by  the  increase  in  damping  ratio  subsequently 
provided.   The  original  system  response  to  a  step  input  is  shown  in  Fig. 22 
with  P  ■  0.147.   Several  other  values  of  beta  were  tested  and  these  re- 
sults appear  in  Fig.  23.   To  verify  both  root-locus  and  parameter  plane 
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predictions  of  instability,  a  test  with  P  =  0.333  was  run  and  the  re- 
sulting response  is  shown  in  Fig.  24.   An  almost  undamped  system  was 
evident,  characteristic  of  a  system  on  the  verge  of  instability  with 
imaginary  roots. 

On  the  basis  of  transient  response  tests,  a  new  operating  point  of 
C  =  2.37,  P=  0.21  was  chosen.   Since  the  thermostatic  level  indicator 
transfer  function  information  was  somewhat  doubtful  as  indicated  in  Ref. 
2,  it  was  felt  that  this  parameter  could  not  be  varied  with  any  meaningful 
results.   The  newly  selected  operating  point  is  marked  with  a  ©  on  the 
parameter  plane,  Fig.  20. 
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IV.   SYSTEM  RESPONSE  TO  VARIATIONS 

Using  the  complete  simulation  program  once  again,  the  system  was 
subjected  to  the  input  shown  in  Fig.  5.   This  time,  the  parameters  used 
on  the  parameter  plane  graphs  of  the  previous  section  were  varied  to 
observe  their  effects  on  the  overall  system.   Note  that  while  a  total 
of  six  parameters  were  investigated,  it  was  recommended  that  only  three 
of  them  be  changed  for  "optimum"  response.   In  the  air  and  oil  loops, 
the  proportional  gain  of  a  controller  was  varied  while  in  the  water  loop, 
one  of  the  combining  relay  gains  was  changed.  For  ease  of  identification, 
these  variables  were  subsequently  called  AIR,  OIL,  and  WAT,  respectively. 

It  was  decided  to  vary  one  parameter  at  a  time  first,  to  observe  the 
effects,  if  any,  of  variable  interaction.  After  this  was  accomplished, 
all  three  were  changed  in  accordance  with  recommendations  made  in  the 
previous  section  and  the  responses  were  again  plotted. 

A.   VARIATIONS  IN  AIR  LOOP 

Table  One  indicates  the  values  of  the  variable  parameter  AIR  for 
the  response  curves  of  Figs.  25  through  29. 


TABLE  ONE 


AIR 

CURVE  # 

1.0  .  . 

1 

2* 

2.0  .  . 

,   3 

.  4 

,   5 

3.5  . 

6 

^Original 

Opera 

ting 

Point 
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Curve  #4  was  the  result  of  using  the  operating  point  recommended 
from  parameter  plane  predictions.   It  was  a  good  deal  faster  than  the 
original  system  as  evidenced  on  the  oil  and  air  flow  responses,  Figs.  26 
and  27.   It  made  little  difference  in  the  two  most  important  variables, 
steam  pressure  and  water  level,  Figs.  28  and  29.   Note  the  near  absence 
of  coupling  of  the  variation  of  AIR  to  water  flow  as  evidenced  in  Fig.  25, 
the  water  flow  response.   The  curves  nearly  coincided  and  were  indistinguish- 
able  as  to  which  curve  was  which. 

B.   VARIATIONS  IN  OIL  LOOP 

Table  Two  indicates  the  value  of  the  variable  parameter  OIL,  for 
the  response  curves  of  Fig.  30  through  34. 


TABLE 

TWO 

OIL 
4.0  .  . 

CURVE  # 

1 
2 
3* 

5.0  .  . 

6.67  .  . 

*0riginal 

Operating 

Point 

As  was  evidenced  from  the  five  response  curves,  Figs.  30  through  34, 
any  proportional  gain  (OIL)  less  than  6.67  slowed  the  system  responses 
considerably.   This  indicated  that  an  increased  gain,  as  predicted  by 
the  parameter  plane  analysis,  was  in  order. 

C.   VARIATIONS  IN  WATER  LOOP 

Table  Three  shows  the  values  of  the  parameter  WAT  for  the  curves 
depicted  in  Figs.  35  through  39. 
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TABLE  THREE 


WAT 
0.21  .  . 

CURVE  # 

1 

2* 
3 
4 

0.147 
0.1  .  . 

0.05  .  . 

*Original 

Ope 

rating 

Point 

All  curves  except  those  in  Fig.  35  practically  coincided  with  one 
another.   Figure  35,  the  water  flow  rate  response  showed  that  the  water 
responded  faster  with  the  new  variable  value  (WAT  =  0.21),  as  predicted 
on  the  parameter  plane. 

D.   VARIATIONS  IN  ALL  THREE  LOOPS 

As  a  final  system  test,  all  three  parameters,  AIR,  OIL,  and  WAT 
were  changed  as  was  recommended  in  the  previous  section  in  accordance 
with  parameter  plane  predictions.  As  a  comparison,  the  old  (original) 
response  was  plotted,  along  with  a  third  curve  representing  the  response 
if  the  parameters  were  moved  in  the  opposite  direction  from  that  pre- 
dicted by  the  parameter  plane.   Table  Four  gives  parameter  values  along 
with  corresponding  curve  numbers  for  Figs.  40  through  44. 


TABLE  FOUR 


AIR 
2.5 
1.67 
1.0 


OIL 
8.0 
6.67 
4.0 


*Orig5nal  Operating  Point 


WAT  CURVE  # 

0.21 1 

0.147 2* 

0.1 3 
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In  each  case,  it  can  be  seen  that  the  response  is  faster  for  curve 
#1  than  for  the  original  system.   This  proves  the  merits  of  design  by  the 
parameter  pl:ne.   The  curve  #3  in  each  case  was  considerably  slower. 

E.    ANOTHER  PARAMETER  VARIATION 

It  had  been  thought  that  the  time  delay  (0.45  sec)  in  the  feed  water 
control  valve  contributed  significantly  to  the  oscillatory  behavior  of  the 
water  flow  loop.   The  system  was  run  using  a  variable  time  delay  in  the 
simulation  to  observe  its  effects.  The  time  delays  used  were  0.45  (original) 
sec,  0.22  sec.  (%) ,  and  0.11  sec.  (%) .  The  response  of  the  water  flow 
loop  is  shown  in  Fig.  45. 

It  can  be  seen  that  the  variation  in  time  delay  did  not  affect  the 
oscillations  but  did  cause  them  to  die  out  more  quickly. 
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V .   CONCLUSIONS  AND  RECOMMENDATIONS 

From  this  work,  it  can  be  seen  that  the  use  of  a  simulation  program, 
particularly  DSL,  was  valuable  in  the  testing  and  analysis  of  a  system  of 
this  type.   The  value  of  other  computerized  analysis  tools  was  also  demon- 
strated.  The  parameter  plane  proved  to  be  most  valuable  and,  even  where 
interaction  was  ignored,  proved  useful  in  the  design  of  boiler  control 
systems. 

The  following  recommendations  are  made  as  a  result  of  this  investi- 
gation: 

1.  Apply  the  simulation  technique  to  other  boiler  control 
systems,  on  other  ships  which  have  the  common  pneumatic  variety  and  on 
the  newer  ships  to  come  which  will  have  a  state-of-the-art  (probably 
fluidic)  system. 

2.  Construct  a  three-dimensional  parameter-space  plot  to 
depict  the  effects  of  varying  the  three  parameters,  OIL,  AIR  and  WAT 
on  the  important  variables,  steam  pressure  and  water  level. 

3.  Investigate  the  effects  of  varying  other  system  para- 
meters or  replacing  certain  system  components,  using  the  original  simu- 
lation. 

4.  Extend  the  original  system  to  be  dynamically  sound  over 

the  entire  operating  range  of  the  boiler  when  other  transient  response  data 
becomes  available. 

5.  Continue  work  on  the  non-interacting  controller  shown  in 
Appendix  B. 

It  is  not  the  belief  of  the  author  that  the  improvements  were  of 
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sufficient  magnitude    to   recommend   changes   to   the  actual   system.      It   is 
hoped,    however,    that   the  methods   and   techniques   used   in   this   thesis  will 
be   of   sufficient  value   to  aid   in   the    further  design  and   improvement   of 
boiler  automatic   combustion  control   systems. 


88 


APPENDIX  A 


BOILER  DESIGN  AND  PERFORMANCE  DATA 


DESIGN 


Steam  Generating  Surface  Area: 
Main  tube  bank 
Screen  tubes 
Front  tubes 
Rear  wall  tubes 
Side  and  roof  wall  tubes 
Economizer 
TOTAL  WATER  SURFACE 
Superheater 

TOTAL  SURFACE  AREA 


Sq.  Ft, 

4975 

933 

40 

158 

274 

7500 

13880 

2390 


16720 


Furance  Volume: 


1270  Cu.  Ft, 


Boiler  Weight: 
Dry 

Water  (Steaming  Level) 
TOTAL  BOILER  WEIGHT 


Lbs. 

230085 

20080 

250165 
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PERFORMANCE  DATA 

Operating  Rate,  7a 

Total  Steam  Generated,  lb/hr 

Drum  Pressure,  psig 

Superheater  Outlet  Pressure , psig 
Temperature,  deg.  F. 

Air  Flow,  lb/hr 

Oil  Flow,  lb/hr 

Anticipated  Efficiency,  % 

Number  of  Burners  in  Operation 


CRUISING 

FULL 

OVERLOAD 

44.7 

100 

120 

117478 

261450 

313800 

1210 

1245 

1265 

1200 

1200 

1200 

950 

970 

965 

145000 

321000 

386000 

8645 

20050 

24100 

87.6 

85.2 

84.2 

4 

6 

6 
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APPENDIX  B 

A  NON-INTERACTING  BOILER  CONTROLLER 

The  utility  of  a  proportional  plus  reset  type  controller  is 
limited  in  an  environment  which  has  many  variables  in  which  inter- 
action is  present.   In  the  CVA-63  boiler  control  system,  for  example, 
a  change  in  steam  flow  affects  the  drum  water  level  as  well  as  the 
pressure.  When  correction  is  applied  to  the  water  flow  rate  to  bring 
the  water  level  back  to  normal,  the  pressure  is  also  affected  which 
causes  an  adjustment  to  be  made  to  the  fuel  flow  rate. 

It  is  desired  to  design  a  controller  which  will  automatically  limit 
the  interaction  between  control  system  variables. 

The  major  input  variable  to  the  boiler  and  control  system  is  the 
steam  flow  rate. 

(la)  APd   =  KPs  AWs   +  Kpf  AWf  +  KplAWl 

(ib)  AX    -  Kxs  AWs    +  Kx-fAWf  +  Kxi  AWi 

These  equations  determine  the  change  in  drum  pressure  (  APd  )  and  water 
level  (  £s.X  )  in  terras  of  steara  flow  rate  change  (AWS),  fuel  flow  rate 
change  (AWf),  ,  ad   feed  water  flow  rate  change  (  AWi.  )  .   The  coeffi- 
cients of  Eqns.  (la)  and  (lb)  are  time  dependent  and  are  determined  from 
boiler  dynamic  analysis. 

The  controller  must  adjust  the  fuel  and  water  flow  rates  in  a  manner 
such  th  .  t  when  the  st^am  flow  rate  changes,  neither  the  drum  water  level 
or  the  pressure  is  changed. 
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The  controller  will  receive  measurements  of  steam  flow  rate  (Ws), 
drum  pressure  (  Pd  )  and  drum  water  level  (X)  from  the  boiler,  compare 
the  latter  two  with  pre-determined  set  points,  Pd  and  X  and  adjust  flow 
of  water  and  fuel  accordingly.   Refer  to  Fig.  46. 

Hence,  the  controller  equations  may  be  written: 

(2a)  AWf    =   M^p(eP)  +  H*x  (ex)-*    Hfs  (AWs) 

(2b)  awl    -  HLp(ep)  +  Hw  (ex)  +  Hcs(aws) 

A  A 

where  ep=Pa-Pa  and  6x-   X-X  are  error  signals.   The  time  dependent  co- 
efficients Kcj  are  known  (Fig.  46)  as  they  represent  boiler  dynamic  be- 
havior.  The  H-coefficients  are  controller  coefficients  and  are  yet  to 
be  determined. 

If  the  equations  (2a)  and  (2b)  are  substituted  into  (la)  and  (lb), 
the  result  is: 

(3a)       APd  =    KPs  AWs  +   Kpf  [  H*p  (ep}  +  H*«  (e*)  +  Hfs  AWs  ] 

+   KpL[Htp(eP)  +  Hu(eO*  H«AWs] 

(3b)       AX=     Kxs  AWs  +    Kxf  LHfP  (eP)  +  H-fx  (e*)  +  Hfs  AWs] 

+    Kxc[HcP(eP)  f- Hc^(ex)+  WcsAWs] 

Rearranging  the  above  equations  and  collecting  terms,  the  above  equa- 
tions can  be  written: 

(^a)     £pd  =  (kPs  +  KpfH*s+  KPtHcs)AWs +  (KPfH+p+  KPLHcP)eP 
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(4b)     z^x  =    (Kxs+KxfHfe+KxiWu)zs^+  (KxfH-fp  4  KxcV4.P)eP 

+  (Kxf  Hfx  +  Kxl  Hu) ex 

Requirements  for  the  non-interacting  controller  force  Afj  and  AX 
to  be  zero  for  any  change  in  GP  or  6<  .   Therefore  the  coefficients 
of  these  terms  must  be  zero  identically. 

(5a)      KPf  HfP  +  Kpi.  Hip  =  O 

(5b)  Kpf  Hfx  *    KP:  Hex    s   0 

(5c)  Kxf  Hfp   +    Kxc  HtP    =  O 

(5d)  Kxf  H-fx    +    KxcHlx    so 

Also  present  is  the  requirement  t  at  when  the  pressure  set  point  is  changed, 
the  pressure  change  must  respond  in  some  predetermined  manner,  say  G(s). 
Then 

(5e)        KPs   +   kf-f  Mrs  +    KPl  Hcs    =   GpCs) 

Similarly,  the  water  level  must  respond  for  a  change  in  water  level  set 
point;  therefore 

(5f)     Kxs  +  Kx-F  rUs  +   Kxc  His  -  GxCS) 

The  six  equations,  (5a)  through  (5f),  contain  six  controller  trans- 
fer functions  (  Htj )  which  are  unknown,  six  boiler  dynamic  transfer 
functions  which  are  known,  and  two  "response  to  set  point  change" 
transfer  functions  which  must  be  specified.   Solving  for  the  six  un- 
knowns yields 
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(6a)  H*s  -    (KvsKpL    -   KPs  Kxc)/  C 

(6b)  Mcs  =   (Kxf  KPs-  KPf  Kx5)/C 

(6c)  H*x=      -(GxKPl)/C 

where  C  =     KPf   Kxd    -    Kpi  Kxf 

(6d)  Hu  =     (GxKpO/C 

(6e)  HfP   -       (GP  Kkc)/C 

(6f)  HcP  -     -(GpKx-O/C 

The  simulation  of  this  controller  was  attempted  using  DSL  in  the 
same  fashion  as  the  original  control  system.  Several  of  the  transfer 
functions  contained  more  zeros  than  poles  and  the  simulation  failed. 

Future  work  on  this  controller  could  involve  the  introduction  of 
excess  poles  in  the  offending  transfer  functions  which  could  be  cal- 
culated to  minimize  any  disturbing  effects  on  the  system. 


95 


APPENDIX  C 

I.    PADE  APPROXIMATIONS 

The  general  form  of  the  fourth-order  Pade  approximation  used  in 
the  simulations  was: 

p-*s  2      Cts)4  -  20(T5)3  4    l80(rs)2-840  rs  +  \(oQ0 


(tsr  +  2o(rsr  +  i8o(rs)2  +840  rs  +  i£8o 

Specific   approximations  were   as    follows    for   the   individual    loops: 
AIR: 

-•05TS„     1.0551  x  lo"5  SA  -     003703  S3  +  .58482  S*  -  47.88S     -t-  IfcSo. 
q  —  — — : 

1.0S5I  <  I0"s    s4  A-   ,  003.701.  S3    +  ,58482  S2   +  47. 88  S    +   1680. 


OIL: 

-  02ss  ,39ot2>cio"feS4  -    3l25xio~3S3+-   .  H25  S^  -  21.  S    +   K80. 

e  " 


.39oGZ  kicks'*  4  /3l2Sxio"3s'i4-    .  1125  S2  +  g|.  s     ■+   1680. 


WATER: 

-.455  04lOOfe2Ss"-    1.8225  S5  -h   3645  S2  -    378.5    +    It80. 

e  ^  _: 

,O4lo0fe25  s4   +■  \.8225S3  +  3t.4S  S2  4-  378.S  +  1480. 


II.   INDIVIDUAL  LOOP  CHARACTERISTIC  POLYNOMIALS: 

The  following  polynomials  were  calculated  using  the  FORMAC  com- 
puter program  and  were  used  in  the  parameter  plane  plotting  routine: 
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AIR:        .l24Gxio      S      +.4&3I  XI0~y  S     +    .$625xlo      S       +  .8598xlo~bS 
+  .4735  x  lo"5s"  +.01329  s'°  +  .  2&I35  S9  4    3.5\34  S8     (  35.736    + 
,52  80xio"7o()  57  4    (258.293  -   .  1691  x  IO~V  4. 52S97 X  Io"7a)  S% 
(|433.6  +  .002364  c(  -   J69lxtcT4A)s5+  (5264.03 -.152950(4. 002364ft)  S^ 
(14502.644   1.637*  -  .  15295  p  )  S3  4    (  12894.84  +  204. 6o3o(  +.  (.637(3)  S2 
(1680.4    2o5fe.32  oc    4  204.603  p)  S    4       2056>,32  A 
OIL:  .2973  X  lo"7  S13       +  .2611    x  I0~4  s'2  +  .0|o3  s"    +    2.144    s'° 

4  197.64   S      +    1765,72   S8  +   I5266.9    S?  +  (  8S5o|.7    4  .1104  X  lo"   cX  )  S& 
4  (233383.7  -  .96  Xlo    ot+.llo4x|o    g\s  +  (286758.5+  .3756c<  _#  %)<  lo    A  )  S4 
4(123122.    -     76.22  c(  4  .3756^)  53  4    (|680.4    6627.8  (X   "76.22  j3)  ^ 

4(29.  856  tX    +66,27.8/3)   S    -4       29.856    8 

WATER:         ,o334fc    S8       4    1.5422    S7    4   32. £4  Sfo    +  (359.64+  .  127  Xio"  j3)  5^ 
4(1924.32   4    ,&38l  x  (o""o(  -.5469  x  lo~3fi)  54    4   (  2748.   -  ,283*  \o~3oi  4 
..0||29|S)s       4-  (2400.72    4  .5672  *   )0~8o(  +    .11718  <3  )  S2     t      (  I68O.  - 
.58824   X  io~To(     +.5208/3)  5     +       .  26l4xiO~    0( 
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TITLE       *****    CVA-63    BOILER    AND    ACC    SYSTEM    SIMULATION    ******* 


*****     INPUT    FORCING    FUNCTION    *******************************   ^ 
AFGEN       STEAM=0.0,Co 0,5,0, 0oG,8o0, 50000. ,llo0,0o0,100o ,0o0      -" 
GS=AFGfcN(STEAM,TIMF)    ' 

*******************  ***************  **********************-*  *** 
****************  A  **  ***********  *********  *  ***************>;**** 

*****    BOILER    PORTION    OF    SIMULATION    ************************* 
PDS=INTGRL(  ICl,-19oE-06*GS) 
*PDF=TRNFR( 1.2 iICA, NUMA,DENA,GF )' 
~PQW=TRNFR(0,3,lCBf  NUMB  ,  DEN  B  ,GW  ) 
PO=PDS+PDF+POW 
PSA=L37E-04*GS 
POUT=PO-PSA 

LF  =  LEDLAG(  I  C 2,  14.2,  lo 42 , 5a 1 1E-06*GF ) 
LF1=INTGRL(1C1,LF> 
LW=INTGRL     UC3f5*02F-07*GW) 
LSA=CMPXPL(  IC4,IC5,2o55,0«63,0.79E-05*GS) 
LSB=INTGRL  (  IC6  *5.0  2E-0  7*GS  >  ' 
LST  =  LSA-LSB"~^^ 
L=LF1+LW+LST 

******  jp.  *********  ******/******************  ***********  ******** 
*****************j>:  ************  **********  ***  *********  ****  A  »  ** 

*****    OIL    FLOW    LOCP    PORTION    OF    SIMULATION '****+************* 
PP=CMPXPL(  IC7,ICR*„042  2,7o09,P0UT*l,>?5)' 
PZP=PRP-PPL * 

PMA=PZP3*OILy  / 

PMB-INTGRL (  IC9,0.14*PZP)  7 
PM=PNA*PMB  ' 

SIGNAL=PQRL-PML    '  y 

PSF=FCNSW(  SIGNAL ♦PQRL,PQRL,PML)/ 


PRF=CMPXPL( IC17, IC18,1.13,lo4B,8*0*PSF) 
PFINT=CMPXPL(IC19t IC20tl.01,2«.485t36»l*PRFLK 
PF=DEL AY ( 10,0.23, PF INT K 
GF=117.*PF- 

************************************************************ 

*******************  **********  ************  ******************* 

*****    AIR     FLOW    LOOP    PORTION    OF    SIMULATION    ****************** 
PEQ=PML-PQ«L^ 

POAA  =  AIR*PEGi/  y 

POAB=INTGPL(IC10,0.024*PEQ)  ' 
POA=POAA+POAB 


PQA=CMPXPL(ATilt_IC12j,0.5,7991  , QA*6o  )  ' 
PQR-TRNFR(1,2  ,  I  CC/NU'MC  ♦  DENC  f  PQAL  >  C 


i/ 


_X,  NU'MC,. 

PQ1  =  TRNF R  (  1  ♦  2  ,  TXT,  NU  NC  ,  DENC  ♦  PO  AL  ) 

GB1INT  =  CMPXPL  (IC13,  I.C14.0.  22  37  ,  5.59,  PQ1* 5 3000.  )' 

-  G  B 1  =  0  E  L  A  Y  { ;L£U-  QJB  i  G  B 1 1  N  T  )  / 

QB1=CMPXPL( IC15,IC16, 1. 534, 0o 42 1 ,0o 444E-03 *GB1 )   »" 

QA=3.0*QB1^ 
**********************************************  ***********^.** 

***.,   ******************************************************** 

*****  WATER  FLOW  LOOP  PORTION  OF  SIMULATION  **************** 
/PSW=.22*PSA-WAT*PWA  x 

PD W=  PR  L  *  PS  L~  P S  W  ^         -t 

GWINT=TRNFR(0,3,[LC_li'NUM0,DEND,P0W)- 
-    GW=DELAY(10,0.457bWINT )   . 
PWA=0. 925E-C4*GW 
PL=2.37*L 

PSL=PS+PL 

**************************************** ******************** 

*********************************************** ******* ****** 
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**  LIMITING  FOR  CERTAIN  PNFU^ATIC  SIGNALS  ** 
-PPL=LIMIT{-20.,?0.,PP} 
PML=LIMIT(  -20.  t20.  ,PM.) 
PORL=LIMIT(-2  0„,20o , PQR  )  ' 
PCAL=LIMIT{-20oCt20*,P0A) 
PRFL=LIMIT(-60o0,60o,PRF)' 
P0AL=LIMIT(-20,»20. ,PQA)  ^ 

************* ******  #*****  *************$.*•*••*;*.}:**************** 
jM*  **************************  :«***t- *************  ************^ 


*****  CONTROL  PORTION  OF  SIMULATION  ************************ 
STOP AG  ICA(2) ,  ICB(3) , ICC(2),  ICO( 3) ,NUMA(2) tDENA(3) , NUMB( i)  ,  - 

DENB<4) ,NUMC(2) »DENC(3) ,NUMD(1 I tDEND(4) 
TABLE  XCA{  1-2)=2^G«C,ICB(  1-3  )  =  3*0„  0,  ICC_LL=21=2*0.»0  ,„  .  , 

ICD(l-3)=3*0o0,... 

NUMA(l-2)=«C62  8f3o 13E-04,DENA( 1-3  )  =71  04  , 1  o  0,  CO,  „  .  o 

NUMB<1I=-4<,42E-06/OENB(  l-4l=10.1t2o54tl.0t0o0t... 

NUMC<l-2)=o064,ln0tDENC( l-3)=.0064f.08,1.0t... 

NUMD(l)=-3. lE04tDEND(l-4)=a816,1.5432,l«204, LO 


C  THESE  PARAMETERS  ARE  FOR  THE 

PAR  AM  0IL=6o67,AIR=1.67fWAT  =  a147 

PARAM  PRP=O.OtPPL=O.OfPS=OoO 

CONTRL   FINTIM=50.,DELT=.05 

INTEG  RKSFX 

END 

STOP 


ORIGINAL  SYSTEM 
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